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ABSTRACT 



Aims. We investigate the contribution of different formation scenarios for type la supernovae in elliptical galaxies. The single degen- 
erate scenario (a white dwarf accreting from a late main sequence or red giant companion) is tested against the double degenerate 
scenario (the spiral-in and merging of two white dwarfs through the emission of gravitational wave radiation). 

Methods. We use a population number synthesis code incorporating the latest physical results in binary evolution, and allowing to 
differentiate between certain physical scenarios (e.g. description of common envelope evolution) and evolutionary parameters (e.g. 
mass transfer efficiency during Roche lobe overflow). The obtained theoretical distributions of the delay times of type la supernovae 
are compared to those which are observed, both in morphological shape and absolute number of events. The critical dependency of 
these distributions on certain parameters is used to constrain the values of the latter. 

Results. We find that the single degenerate scenario alone cannot explain the morphological shape of the observational delay time 
distribution, while the double degenerate scenario (or a combination of both) can. Most of these double degenerate type la supernovae 
are created through a normal quasi-conservative Roche lobe overflow followed by a common envelope phase, not through two suc- 
cessive common envelope phases. This may cast doubt on the use in other studies of analytical formalisms to determine delay times. 
In terms of absolute number, theoretical supernova la rates in old elliptical galaxies lie a factor of at least three below the observed 
ones. We propose a solution involving the effect of rotation on the evolution of intermediate mass binaries. 

Key words, supernovae: general - binaries: close - stars: white dwarfs - galaxies: elliptical 
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1. Introduction 

Type la supernovae (SNe la) are among the most powerful events 
observed in the universe. Their seemingly identical nature has 
led them to be most useful distance indicators, or standard can- 
dles, in cosmology, while their chemical enrichment of the inter- 
stellar medium is of paramount importance for the chemical evo- 
lution of galaxies. Large amounts of the chemical element iron 
(Fe) are released during a SN la. As an often quoted example, 
the present mass fraction of this element in the solar neighbor- 
hood could not be explained without the inclusion of SNe la in 
chemical evolution models. 

The current consensus is that SNe la originate from the ther- 
monuclear disruption of a white dwarf (WD) which is no longer 
able to support its own mass through degenerate electron pres- 
sure, as a resul t of approa ching or exceeding the Chandrasekhar 
limit (see e.g. iLiviol [200lb . However, many open questions re- 
main: whether only carbon-oxygen (C-O) WDs are involved, 
whether the disruption is a detonation or a deflagration, whether 
the event takes place at the Chan drasekhar mass or at a sub- 
Chandrasekhar mass, etc. (see e.g. lBranch et alJ ll995). Perhaps 
most importantly, even the exact progenitor system of a SN la 
has not been unambiguously established. There exist different 
scenarios, some of them involving one degenerate object accret- 
ing from a less evolved companion in a multiple star system, 
others assuming that SNe la are caused by the merger of two 
WDs. All of these elements can have important consequences 



for the applications of SNe la, especially for their use as stan- 
dard candles. 



In this work, we differentiate between the two most in- 
voked formation scenarios for SNe la: the single degenerate 
(SD) and double d egenerate (DD) model. In the SD model (see 
e.g. lNomotol [T982) a WD accretes matter through mass transfer 
from a companion in a close binary system. The companion can 
be either a late main seq uence (MS) star or a red giant (RG). In 
the DD model (see e.g. Ilben & Tutukovlll984: IWebbinki ri984) 
the SN la is caused by the merger of two WDs, the combined 
mass of which equals or exceeds the Chandrasekhar mass. Their 
merger is the result of a spiral-in caused by the emission of grav- 
itational wave radiation (GWR). The validity of this scenario has 
often been questioned, since the exact detonation mechanism 
is unclear. It has been argued that a merger may instead result 
in an off-center carbon ignition and the formati on of a oxygen- 
neon- magnesium WD or a neutron star (s ee e.g.LS aio & Nomoto 
1998). However, by including rotation iPiersanti et al.l ([2003) 
have shown that the accretor star can cool down, causing a 
greater expansion of its outer layers and a lowering of the critical 
angular velocity. This limits the mass accretion rate and prevents 
the occurrence of a gravitational instability, allowing mass accu- 
mulation up to the explosive central ignition of carbon, leading 
toaSNIa. 
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1.1. The delay time distribution 

The question that will be addressed in this work is which of these 
two scenarios is the main contributor to the SN la explosions ob- 
served in nature. In order to do this, a comparison with observa- 
tions is obviously necessary. Observational studies of the distri- 
bution of delay times of SNe la can be used to constrain fheo- 
retica l models for their progenitors and formati on scenarios (see 
e.g. iMannucci et al.ll2006h . iTotani etal] d2008l) selected 65 su- 
pernova candidates at redshifts from 0.4 to 1 .2 (corresponding to 
a light travel time of 0. 1 to 8.0 Gyr), of which they demonstrate 
that at least 80% must be SNe la. These events were selected 
from old galaxies of the elliptical type, in which star formation 
has long ceased almost completely. Hence, the selected samples 
are equivalent to (passively evolving) starburst galaxies for the 
determination of the delay time of SNe la, and can be used for 
that purpose. The thus obtained observational delay time distri- 
bution (DTD) can then be directly com pared to the theore tically 
predicted DTDs for starbursts. While Totani et al. (2008) make 
the implicit assumption that the DTDs are independent of metal- 
licity, the samples are expected to be contained within the range 
1.0 to 2.5 times solar metallicity, so any metallicity effect should 
be limited and certainly less than an order of magnitude (as will 
be shown later). Also, it is assumed that the average age of local 
elliptical galaxies is 1 1 Gyr, and for tho se the authors adopt th e 
SN la rate observed for such galaxies bv lMannucci et al. (2005). 
A striking feature of the thus obtained observational DTD is that 
it follows a power law: it decreases inversely proportional to time 
(galaxy age). 

The only caveat is that the observations are presented per unit 
delay time per century for a single starburst population whose 
total K-band luminosity L^o at an age of 11 Gyr is 10 10 Lk,©, 
while all population synthesis models give SN la rates in num- 
ber of events per unit time for a single starburst with a given 
total initial mass. This implies that a conversion is required be- 
tween SN la rate p er K-band luminosit y (SNuK) and SN la rate 
per mass (SNuM). [Totani et al. (2008) adopt a conversion fac- 
tor obtained from a spectral energy distribution template for an 
exponentially decaying star formation rate with a timescale of 
0. 1 Gyr. This factor depends on the adopted metallicity Z and 
initial mass function (IMF), and values are given for different 
choice s. The cases assuming solar metallicity and the Chabrier 
(2003 1) IMF best match our use (see Sect. 2) of Z = 0.02 and 
the iRroupa et alj 0993) IMF. Following personal communica- 
tion with Totani, we take a value based on a combination of these 
models: 



M, 



= 1.56 



L 



K.O 



(1) 



In this equation, M, is defined as the integrated star forma- 
tion rate up to the age of the galaxy, which in the case of a single 
starburst is nothing else than the total initial mass. Since this con- 
version factor is taken to be constant for the entire redshift range, 
any change (or uncertainty) in it will only influence the absolute 
values of the SN la rate in SNuM, and not the morphological 
shape of the DTD in those units. 

The question addressed in this work will thus pertain to the 
dominant formation scenario in the case of elliptical galaxies. As 
a novelty, we will specifically investigate the influence of the de- 
scription of mass transfer efficiency during Roche lobe overflow 
(RLOF). In Sect. 2, we give a brief description of our population 
number synthesis code, and we elaborate on the new physical 
effects which have been incorporated as well as on the details 
of systems leading to a SN la. Sect. 3 contains our new results 



on the progenitors of SNe la, as well as a short elaboration on a 
special kind of merger. In Sect. 4, a comparison is made between 
the results of our code and those published by groups undertak- 
ing similar studies. 



2. Population number synthesis code 

For our study, we use the Brussels population number syn- 
thesis code developed over the years at the Astrophysical 
Institute of the Vrije Universiteit Brussel. A very extended 
review of the elements contain ed in this code is given by 
iDe Ponder & Vanbeverenl (l2004). We repeat here only the most 
important points and the changes which have been implemented 
since then, reflecting recent developments in the description of 
certain physical processes. 

Our code starts from an instantaneous starburst normalized 
to 10 6 M with Z = 0.02 and a binary fraction of 1 00%. By 
default the stars are drawn from alKroupa et al.l d 19931) IMF, flat 
mass ratio distribution and Abt (1983) initial orbital separation 
distribution. The implications of other choices will be discussed 
in Sect. 3. Supernova rates are given in the output as number of 
events per Gyr as a function of time elapsed since starburst. 



2. 1. Treatment of mass transfer 

As stated above, a main new factor in this investigation is the 
influence of the description of the RLOF. Whether or not this 
process is conservative is indicated by the parameter ft. 

In the conservative case, (3 — 1, the total mass and angular 
momentum of the system are conserved during RLOF. Matter 
is transferred from the initially more massive star (Mi) to the 
initially less massive star {Mi), whereby the latter accretes all 
the matter lost by the donor. This process can continue after the 
donor star has become less massive than the accretor star, a phase 
known as mass ratio reversal. Angular momentum is exchanged 
between the two components, and with the orbit. The resulting 
orbital period is given by 



Pi 
Pi 



M 



iM 2i \ 3 



M lf M 2f / 



(2) 



The indices in this formula, i for initial and f for final, are 
with respect to the RLOF phase. This ratio is smaller or larger 
than unity depending on whether mass transfer continued until 
after mass ratio reversal. 

When mass transfer is non-conservative, the following rela- 
tion holds: 



dM 2 /dt 



dMi/df 



< 1. 



(3) 



/3 is thus the fraction of material lost by the donor which is 
accepted by the accretor. The amount of mass lost from the sys- 
tem is readily found by integration over the time during which 
mass loss takes place. To determine the amount of angular mo- 
mentum lost, and thus the final orbital period, however, requires 
knowledge of the way in which the non-accreted material leaves 
the system. The two most common approaches, which will be 
considered in this work, is that either the matter is lost with 
the specific angular momentum of the accretor, or with the or- 
bital angular momentum of a corotating point at the second 
Lagrangian point (L2). The first assumption is justified if the 
mass loss can be treated as an enhanced stellar wind, or at least as 
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a process which is symmetrical in the equatorial plane of the ac- 
cretor. There is however no physical mechanism which is known 
to make this possible in normal intermediate mass binaries. 

The formulas thus obtained for orbital period variation 
during non-conservative RLOF are the following (see e.g. 
IDe Ponder & Vanbeverenl2004l) . 

ForO <p< 1: 



P, \M li + M 2i \Mu \M 2i 



-3[V^ + l] 



(4) 



P 



For 13 = 0: 
/ Mu + M 2{ 



Mu + M 2i \Mu 



(5) 



In the case where mass loss through the L 2 point is assumed, 
the parameter r\ is determined by the distance from this point to 
the center of ma ss, expressed in terms of the orb ital separation. 
It was shown bv lDe Ponder & VanbevererJ d2004l) that a typical 
value is r\ — 2.3. When the angular momentum loss is taken 
to be that specific to the accretor, r\ can be much smaller than 
unity. The standard model will assume the former value of 2.3, 
the implications of other choices being briefly discussed. 

When mass transfer becomes dynamically unstable, the bi- 
nary will be enveloped in the extended atmosphere of the ini- 
tially most massive component, resulting in what is known as a 
common envelope (CE) phase. The friction due to viscosity of 
the rotating pair within the envelope causes a significant loss of 
rotational energy, which leads to a strong reduction of the or- 
bital period. This process is known as spiral-in. If the period is 
sufficiently reduced, both stars will merge. Otherwise, if there 
is sufficient binding energy released by the inspiraling pair, the 
envelope may eventually be ejected, again yielding a separated 
binary. Due to the short timescale of this process, it is assumed 
that during a CE phase the accretor does not gain any appreciable 
amount of mass, i.e. /3 — 0. 

Different models exist for the description of en ergy conver- 
sion during CE. The a-scenario bv lWebbinkl (fl984) depends on 
the efficiency a with which orbital energy lost by the inspiraling 
stars can be converted into kinetic energy for the expulsion of 
the envelope. The change in orbital separation is given by 



M u (M u - M U ) 



AR 



Roche 



MuM 2i M n M 2i 



2A f 



2Ai 



(6) 



In this equation, PRoche is the Roche radius of M\, while A 
is characterizing the density structure of this star's extended en- 
velope and is dependent on whether this envelope is radiative or 
convective. By definition, a must lie between and 1, while it 
is also conceivable that the value is dependent on the nature of 
the object which is performing the spiral-in. The resulting or- 
bital period variation is obtained using Kepler's third law. When 
we apply this scenario in our code, the standard choice will be 
that aA = 1, wh ere the Roche lobe is calculated using the fit by 
Eggl etonl d!983h . The implications of different choices will be 
described below. 

An alternative descriptio n of the CE evolution phase is given 
bv lNelemans & Toul d2005l) . The change in angular momentum 
is described by 



A.J 



J 7 Aftotal 



AMtotai Mu - M u 

= r- 



M Vl + M 2i 



(7) 



where y is a parameter with a typical value around 1 .5. Using 
this formalism, the change in orbital separation is as follows: 



Af 
A] 



~Mif 



Mu + M 2i 
Mu + M 2i 
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Mu - Mu 
Mu + M 2i 



(8) 



This equation can again be readily converted into one for 
the evolution of the orbital period using Kepler's third law. 
The authors argue that this formalism is to be preferred in 
systems which will lead to the formation of a double WD 
or WD+MS/RG binary. In brief, this y-scenario starts from a 
conservation of energy and a balance of angular momentum, 
whereas the a-formalism assumes a conservation of angular mo- 
mentum to arrive at a balance of energy. 

2.2. Single degenerate scenario progenitors 

Hachisu et al ] (l!996h introduce a SD SN la progenitor channel 
consisting of a WD accreting hydrogen rich material from a low 
mass RG. High mass accretion rates lead to the development 
of a strong wind, which stabilizes the mass transfer rate. The 
WD thus steadily accretes matter, which is first burned into he- 
lium, without the fo rmation of a CE, allowing i t to reach the 
Chandrasekhar mass. lLi & van den Heuvell(ll997l) e xtended this 
scenario to MS companions. Hachisu et al. (1999b) give a sce- 
nario which takes into account both possibilities for Z = 0.02. 
For different WD masses, regions in the orbital period - com- 
panion mass parameter space are presented which will eventu- 
ally result in a SD SN la. There are two such 'islands' in the 
(P, M compan j on )-plane: one for (late) MS companions, ranging 
from about 2 to 3 M Q companions with orbital periods of about 
0.4 to 6 days, and a second one for RG companions, ranging 
from about 1 to 3 M Q and 30 to 300 days. The first channel 
is known as the WD+MS channel, the second as the WD+RG 
channel. If at any time between the formation of the first and 
second WD in our potential progenitors, the system enters one 
of these contours, it is assumed that it will eventually result in a 
SN la through the SD scenario. The explosion is taken to occur 
at the time that the mass transfer toward the WD starts, as the 
extra time which is required in r eality is neglig i ble co mpared to 
the nuclear lifetime of both stars. Hachisu et al. 1 (120081) introduce 
the influence of the mass stripping effect on the WD+MS chan- 
nel. This effect, resulting from the interaction of the WD wind 
with the surface of the companion, attenuates the mass transfer 
rate from the MS star onto the WD. MS companions can thus 
be more massive than previously expected, and still avoid a CE 
phase, allowing them to remain a potential SD SN la progenitor. 
The fraction of mass stripped away by this wind is 



dM, 



dr 



Mnp dM w ; n d 
= Cl- 



(9) 



ci is a parameter depending on, among others, the wind ve- 
locity, the efficiency of conversion from kinetic energy to ther- 
mal energy, and the geometrical factor of the surface hit by the 
wind. The value of this parameter is quite uncertain, and best 
fit models based on observations yield a possible range from 1 
to 10. A brief discussion on the implications of varying this pa- 
rameter will be given, while a moderate value of c\ = 3 will be 
adopted for the standard parameter set. In that case, the WD+MS 
channel contour typically extends up to a companion mass of 6 
M Q , instead of 2 Mq witho ut the effect. 

iKobavashi et al.l d 19981) investigated the influence of metal- 
licity on the SD scenario. They find that if the accretor's iron 
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abundance is too low, its wind is too weak to allow stable mass 
accumulation up to the Chandrasekhar mass, and hence no SN 
la can occur. The area in the parameter space which leads to 
SD events thus decreases with metallicity. We use their progen- 
itor regions i n the (P, M comDan j on )-plane fo r Z = 0.004 together 
with those of lHachisu et al.l d!999bl 120081) to interpolate for any 
metallicity. 



2.3. Double degenerate scenario progenitors 

In the DD scenario, it is assumed that every WD merger ex- 
ceeding 1 .4 M results in a SN la. In our code, which follows 
the evolution of the progenitor in detail throughout its existence, 
without the use of an analytical formalism, there are two typi- 
cal formation channels which lead to a SN la. We have termed 
them the RLOF channel and the CE channel. A description as 
well as a typical example of both is given below, since experi- 
ence has shown that there exists a great diversity in evolutions 
obtained with different codes (see also Sect. 4). The numbers 
in this subsection are given for a somewhat stricter subset of the 
progenitors of DD SNe la, i.e. those in which both merging WDs 
are of the C-O type. If this requirement is waived, the parame- 
ter space of progenitors is slightly extended. The implications 
of this limitation on the eventual DTDs will be discussed in the 
Results section. 



2.3.1. The Roche lobe overflow channel 

In this channel, the progenitor system at birth consists of a pri- 
mary typically between 3.7 and 7.1 M , a mass ratio between 
0.35 and 1.0 (high mass ratios being strongly favored in case of 
all but the highest primary masses), and an orbital period of up to 
100 days. After the most massive component evolves off the MS 
and traverses the Hertzsprung gap, a first phase of mass trans- 
fer will take place. The envelope of the donor is radiative at that 
time, resulting in a dynamically stable mass transfer on a nuclear 
or thermal timescale, which is a canonical RLOF. The standard 
model is to treat this RLOF conservatively (J3 = 1), i.e. all matter 
lost by the donor is accreted by the other star. The mass transfer 
will continue until well after mass ratio reversal, resulting in a 
system with a more extreme mass ratio (where the donor is now 
the least massive star) and an increased orbital period, typically 
a few hundred days. The evolution of the accretor under the ef- 
fect of ma ss accretion (i.e. re j uvena tion) has been explicitly cal- 
culated by IVanbeveren et alJ (Il998l). This ha s been done with a 
standard accretion model (see Neo et alll 19771) if matter directly 
impacts the surface of the star, a nd under the assumption of ac- 
cretion induced full mixing (see Vanb everen & De Loordfl9"94T) 
if a Keplerian accretion disk is formed. During the collapse of 
the former donor to a WD, a fraction of its mass will be expelled, 
resulting in a somewhat larger period still. After the initially least 
massive star eventually also leaves the MS (a timescale which 
is critically affected by the amount of mass it received during 
the first mass transfer phase), it will also increase in radius until 
it exceeds its Roche lobe. Thus, a second mass transfer phase, 
toward the WD, will be initiated. Due to the extreme mass ra- 
tio of the system and the fact that the accretor is a WD, this 
mass transfer will become dynamically unstable, resulting in a 
CE phase leading to a spiral-in. Eventually, the envelope will be 
ejected and a double WD binary will emerge. This system con- 
sists of two WDs with a mass around one solar mass each and an 
orbital period which has decreased dramatically to a few hours 
during the spiral-in. Constant emission of GWR will cause the 



two components to slowly spiral in over the course of a few bil- 
lion years. Eventually, they will merge, and since they together 
exceed the Chandrasekhar mass, will under our assumptions re- 
sult in a SN la. The typical timescale of this event is the sum of 
the formation timescale of the double WD binary and the GWR 
spiral-in. This is typically on the order of a few Gyr, but can be 
lower for systems which had very short initial periods, resulting 
in a delay time varying between 0.22 Gyr up to the Hubble time 
and beyond. 

If the same initial system is considered, but the first mass 
transfer phase is taken to be totally non-conservative (J3 — 0) the 
system will merge before it can result in a double WD binary, 
and thus there will not be a SN la explosion. 



Example We start from a system on the ZAMS consisting of 
4.0+3.6 M stars with an initial orbital period of 5.0 days. After 
the first mass transfer phase at 0.20 Gyr, which is a canonical 
RLOF and is considered to be conservative, a 0.65+7.0 M Q sys- 
tem emerges with a period of 160 days. Just prior to the origi- 
nally least massive star filling its Roche lobe, the other star has 
become a WD with a mass of 0.57 M and the period has in- 
creased to 170 days. This second mass transfer phase at 0.26 Gyr 
results in a CE evolution (described with the a-formalism), even- 
tually yielding a system of 0.57 + 0.88 M WDs with a period 
of 0.19 days. This requires 3.1 Gyr of GWR in order to merge, 
resulting in a DD SN la after 3.4 Gyr in total. The evolution of 
this system is represented in the left panel of Figure 1 . Using the 
same initial conditions but assuming a totally non-conservative 
RLOF as first mass transfer phase, the system already merges 
during that phase, and thus does not result in a SN la. 



2.3.2. Common envelope channel 

Alternatively, one can start with a system which is quite similar 
in initial masses, i.e. a primary mass between 5.3 and 7.9 M and 
a mass ratio between 0.84 and 1.0, but with a much larger initial 
orbital period, starting at 200 days. The large orbital separation 
will ensure that the most massive star does not exceed its Roche 
lobe until after it has developed a deep convective envelope. The 
result is that this mass transfer phase will be dynamically un- 
stable, and will already here result in a CE spiral-in. This will 
decrease the orbital period by about two orders of magnitude. 
After the envelope is ejected, a binary with a less extreme mass 
ratio than in the RLOF scenario (the accretor has not gained any 
material due to the fast timescale of the mass transfer) and an or- 
bital period of only a few days is obtained. When the originally 
most massive star has become a WD and the other has started 
evolving off the MS, a second mass transfer phase will be initi- 
ated. Just as before, this will result in a CE phase and spiral-in. 
Systems which do not merge during this phase will eventually 
emerge as a double WD binary of the order of one solar mass 
each, with an orbital period of only a few tens of seconds. GWR 
has to work only a few tens of thousands of years in order to 
let such a system merge, which means that the resulting SN la 
is not delayed much more than the formation timescale of the 
double WD binary. This typically amounts to somewhat more 
than a hundred million years, but can vary from 55 to 180 Myr 
depending on the initial conditions. 



Example An initial system with 6.0+5.4 M stars and an initial 
orbital period of 1000 days encounters a CE phase during the 
first mass transfer at 83 Myr, emerging as a 0.91+5.4 system 
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RLOF with ft = 1 (SN la at 3.4 Gyr) 



4.0+3.6 M @ 
P= 5.0 d 
t= 




0.57+6.7 M 
P = 170 d 
t = 0.26 Gyr 



o. 



0.65+7.0 M @ 
P = 160 d 
f = 0.20 Gyr 



o 



0.57+0.88 M @ 
P= 0.19 d 
t = 0.26 Gyr 



CE (SN la at 97 Myr) 



6.0+5.4 M 
P= 1000 d 

t = 



0.80+5.4 M s 
/>=5.6d 
t = 96 Myr 



0.91+5.4 M s 
P= 5.6 d 
£ = 83 Myr 



0.80+0.76 M @ 
P = 700 s 
£ = 96 Myr 



Fig. 1. A graphical representation (not to scale) of the evolution of a system leading to a DD SN la through the RLOF channel (left 
panel) and the CE channel (right panel). 



with a period of 5.6 days. The second mass transfer phase at 
96 Myr, toward a 0.80 M WD, also results in spiral-in. The 
eventual result is a 0.80+0.76 M WD binary with a period of 
700 s, which needs 650 kyr of GWR and results in a DD SN 
la 97 Myr after star formation. This evolution is rendered in the 
right panel of Figure 1 . 

2.4. White dwarf plus main sequence or red giant mergers 

If at any time during the evolution, after the originally most 
massive star has become a WD but while the other is still non- 
degene rate, the components me rge, the result is a priori un- 
known. Sparks & StecheJ (1 19741) suggested that this could lead 
to a supernova explosion. Our code calculates the number of 
such mergers for which the combined mass of the WD and the 
companion's core exceeds the Chandrasekhar mass. As the MS 
or RG star still has hydrogen in its outer layers, we will tenta- 
tively classify such events as type II supernovae. It is however 
not yet known whether this is the case in nature. 

3. Results and discussion 

3.1. Single versus double degenerate scenario 

The DTDs which were obtained with the described code are ren- 
dered in Figure 2. This shows the DTDs using the a-scenario for 
CE evolution (a - I) and with fi = 1, for both the SD and DD 
model. Also re presented are the ob servational data points, with 
error bars, from lTotani et alJ (120081) . The observations have been 
converted from SN la rate per K-band luminosity into SN la rate 
per ma ss using Eq. (1). This includes the one by Mannucci et al. 
(2005) for local ellipticals, which are assumed to have an age of 
11 Gyr. 

The most important observations are the following: 
• As is clear from the figure, during most of the time after star 
formation, the DD DTD is dominant over the SD one, at later 
times by an order of magnitude or more. For much of the first 
1 .4 Gyr, the SD DTD shows a spike which causes it to be domi- 
nant during this time, also by up to about an order of magnitude. 
After these early events however, the SD DTD decreases rapidly, 
then levels off and after about 8 Gyr disappears completely. This 
cutoff time is determined by the nuclear lifetime of the least mas- 
sive companion stars that can lead to a SD SN la. The SD DTD 



alone is thus clearly incompatible with the morphological shape 
dictated by the observational data points. 

• The DD DTD follows a decay inversely proportional to time, 
the shape of which is in fair agreement with that of the data 
points. A DTD made up of the sum of events caused through 
both the SD and DD scenario matches most observational points 
even better in morphological shape, apart from those at very 
early times. It is indeed possible that both models should be con- 
sidered together, since they originate from two distinct groups of 
progenitor systems, whose initial parameters do not overlap. 

• In terms of absolute number of events, neither the SD sce- 
nario, DD scenario, or a combination of both is able to produce 
enough SNe la to match the observed values, especially at later 
times. At the 1 1 Gyr point the discrepancy amounts to a factor of 
three, a number which would even increase if the initial binary 
fraction is taken smaller than 100%. Similar conclusions have al- 
ready been reached by o ther groups perform ing theoretical DTD 
studies (see Sect. 4), e.g. lRuiter et a l. (2009) who found a rate al- 
most t en times smaller than the SNuM given by Mannucc i et alJ 
(2005). A possible explanation is that the constant conversion of 
the observational results from SNuK into SNuM may introduce 
large systematic errors, which is visible from the uncertainties 
on the c onversion factor as wel l as t he difference in its values 
given by iMannucci et alJ d2005h and iTotani et all d2008h . As a 
result, the remainder of this section will mainly concentrate on 
a comparison between the morphological shapes of theoretical 
and observational DTDs, not on their absolute values. 

• We find that nearly all SD events are created through the 
WD+MS channel described above. 

• A great majority of DD SNe la (about 80%) are formed 
through a RLOF phase followed by a CE phase, as opposed to 
two successive CE phases. This will be further discussed in the 
next subsection. 



3.2. Mass transfer efficiency 

In Figure 3, the influence of mass transfer efficiency is inves- 
tigated. It shows the DD DTD for different values of accreted 
mass fraction /3. The figure shows that in the case of totally non- 
conservative RLOF (J3 = 0), the DTD drops dramatically after 
0.2 Gyr, leaving virtually no events beyond 0.5 Gyr. This directly 
implies that in the case of a non-conservative first mass transfer 
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Fig. 2. DTDs of a 10 6 M starburst obtaine d for the DP (solid ) and SD (dotted) scenario w ith (3 
CE-evolution. Observational data points of Tota ni et al.l d2008l) and lMannucci et al] d2005h . 
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phase, there will be no DD SNe la with a sizable delay time. 
Since non-conservative RLOF excludes systems going through 
such a phase as SN la progenitor, the events which do remain in 
this case (all with short delay time) must be produced through 
the double CE channel. On the other hand, the events which are 
absent in the case of (3 — 0, but account for about 80% of all 
DD SNe la in the case of conservative RLOF, are those created 
through a RLOF phase followed by a CE evolution. These results 
can be understood in terms of the two typical evolution channels 
leading to a DD SN la discussed before. It was shown that the 
RLOF channel (one RLOF followed by one CE phase) produces 
events with a delay time of multiple Gyr, but is unable to pro- 
duce any SN la if RLOF is non-conservative. The CE channel 
(two successive CE phases) on the other hand is independent of 
P, but can only produce events during the first few hundreds of 
millions of years after starburst. 

In order for the model to be consistent in DTD shape with 
observations, a nearly conservative mass transfer is necessary, 
which means that j3 needs to be close to one overall. Further in- 
vestigation has pointed out that, on average, f3 > 0.9 is required 
for the theoretical model to remain marginally consistent in DTD 
shape with the observations under the current assumptions. This 
is also visible from Figure 3. An important consequence is that 
taking the formation timescale of a double WD binary equal to 
the MS lifetime of the least massive component is unjustified. 
This practice, commonly adopted in studies which make use of 
analytical formalisms instead of detailed evolution, is namely 
based on the unfounded assumption that the accretor is not af- 
fected by the mass transfer process, either because it is totally 
non-conservative or a CE evolution (which is by definition taken 
to be too fast for the accretor to gain any mass). 

The description of angular momentum loss does not affect 
this quasi-conservative model much, since there is only very lit- 
tle mass which leaves the system. However, when calculations 
are made with lower values of (3, the DTD is critically affected 
by the choice for this description. Mass loss through the corotat- 
ing ring causes a DTD which descends much steeper than with 
the assumption of specific accretor angular momentum loss. This 
implies that the latter description allows a broader range of (3 
to be compatible in DTD shape with the observations (see also 
Sect. 4). 



The SD DTD has also been calculated with (3 = 0, but in that 
case there is little difference with the result shown in Figure 2: 
the DTD drops only very slightly, since the value of /3 is of only 
small influence on the number of progenitor syste ms which will 
at som e point lie within the zones identified by lHachisu et al.l 
(l2008h . This is because some 90% of these systems have gone 
through a CE phase during the first mass transfer, and are thus 
not affected by the degree of conservatism of stable RLOF. 

All these results were obtained under the assumption that all 
WD mergers exceeding 1 .4 M result in a SN la. If this assump- 
tion is changed into including only the merger of C-O WDs sat- 
isfying the same criterion, the absolute number of DD SNe la de- 
creases by about 40%, and the separation in time between RLOF 
and CE events thus becomes more distinct. However, none of the 
above conclusions are altered in any way, since the overall shape 
of the DTDs does not change. 

3.3. Influence of other parameters 

The influence of the c\ parameter in the mass stripping effect is 
investigate d next. So far, c\ — 3 represented the standard model 
favored bv lHachisu et al.l d2008l) . When c\ is taken equal to 0, the 
effect is essentially turned off, and the earliest and highest spike 
in the SD DTD disappears. This makes the absolute number of 
SD SNe la during early times (< 0.5 Gyr) no longer exceeding 
those of the DD channel, effectively making the latter dominant 
over the entire time range. Taking c\ > 3, up to 10 according to 
the empirically obtained values, further enhances the peak which 
is present for c\ = 3, but not by a significant amount. 

Results discussed so far have been obtained with the a- 
scenario of Webbink (1984) for CE evolution. When the y- 
scenario by Nelemans & Toutl (120051) is used (y = 1 -5), the re- 
sults obtained for f3 — 1 are shown in Figure 4. It is obvious that 
the SD DTD still has a morphological shape which is incom- 
patible with the observational data points: it drops away much 
too soon and too fast. Whereas the DD DTD lies almost an or- 
der of magnitude lower in absolute rate than when using the a- 
scenario, this result can by no means be rejected based on a com- 
parison with the morphological shape of the observational DTD. 
An interesting point is that when f3 is taken smaller than one, 
the number of SNe la in both channels decreases (as with the 
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a-scenario), albeit the DD DTD does not disappear as dramati- 
cally as before. For /? = it drops by only about 10%, rendering 
a broad range of /? compatible in DTD shape with observations. 
The reason is that with the y-scenario, there are much more sys- 
tems which are formed through a double CE phase, and are thus 
not affected by the value of p. It is very important however to re- 
mark that the validity of the y-scenario is questionable due to the 
fact that almost no events predicted by it result from the merger 
of two C-O WDs. If progenitors are restricted to such systems, 
adoption of the y-scenario will thus hardly yield any DD SNe la. 

When the c-scenario is considered again, but this time with 
a - 0.5 instead of a — 1, the shape of both DTDs and the ab- 
solute number of events change, but not in a way to alter con- 
clusions. The DD DTD comes to lie higher (+60% in absolute 
number of SNe la), but is still consistent in shape with observa- 
tions, possibly in combination with the SD DTD. The latter one 
drops slightly compared with a — 1 (-30%), but does not change 
in overall shape. The dependence of the DTDs on the value of (3 
is similar to that in the case of a = 1, but with an even earlier 
cutoff for the DD DTD in the case of p = 0. 

Also the influence of initial mass ratio distribution is investi- 
gated. Using a different distribut ion than the flat one use d so far 
(e.g. the ones bv lGarmanv et alj|1980t iHogeveenll 19921) results 
in a change of the absolute number of SNe la in both scenarios, 
and can also slightly affect the morphological shape of the the- 
oretical DTDs. However, the change is too insignificant to alter 
the conclusions above. The same is tr ue for the choice of IMF: 
until now, we have used the one bv iKroupa et al.l dl993h . and 
taking another one (e.g. following ISalpetei 1 19551) does have an 
effect on the absolute number of events, but not on the shape of 
the DTDs or the conclusions made so far. 

A final parameter which is scrutinized is the metallicity of 
the galaxy. So far, calculations have been made with a solar Z = 
0.02, since this is the order of magnitude at which the available 
observations have been made. If we take a Magellanic Cloud 
metallicity of Z = 0.002, the DD DTD is not much affected in 
shape and increases by about 50% in number, whereas the SD 
DTD significantly drops over its entire range, leaving only 3% 
of the absolute number of SD SNe la. This is a direct result of the 
fact that the areas in the (P, M compan i „)-plane of systems which 
will result in a SD SN la decrease dramatically with metallicity. 



An interesting consequence is that very old galaxies with Z close 
to zero are unable to produce a sizable number of SNe la through 
the SD scenario. 

In this respect, it should be remarked that our previous con- 
clusions are not compromised by the combination of this metal- 
licity effect and the fact that the observational data points were 
taken from samples with a metallicity of up to 2.5 times the solar 
value. The reason is that for Z = 0.05 the theoretical SD DTD 
stays within an order of magnitude of the Z = 0.02 one, and its 
shape also strongly deviates from that of the observational data 
points. 



3.4. White dwarf plus main sequence or red giant mergers 

Figure 5 shows the DTD for the merger of a WD and a MS star 
or RG. This demonstrates that for both extreme values of f3, these 
(potential) special type II supernovae start very soon after star- 
burst (about 20 Myr) and extend over a few hundred Myr. They 
originate from systems with a wide range of initial primary mass 
(above 4.5 M ), mass ratio (between 0.2 and 1.0) and orbital pe- 
riod. In the case of conservative RLOF (P — 1), about one quar- 
ter of the events are formed through two successive CE phases. 
These are the only ones which remain for p — 0. 



4. Comparison with previous studies 

The results of our population number synthesis can, after a sim- 
ple conversion of units, be directly compared to results obtained 
by other groups in similar theoretical population studies. As 
sketched above, population synthesis models are very sensitive 
to a large number of input parameters, such as mass transfer ef- 
ficiency, CE description, angular momentum loss mechanism, 
initial mass, mass ratio and orbital period distribution, and oth- 
ers. In order to carry out the following comparisons, we have 
equipped our code with the assumptions, parameters and initial 
conditions of the other groups, as far as the code technically al- 
lows. The descriptions of the other codes are given in the next 
subsections, while the comparison with ours is made in Sect. 4.5. 
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4.1. Yungelson & Livio 

The first study with whi ch a comparison is made is that of 
Yung elson & Liviol d2000l) . The authors give, among others, a 
study of the DD scenario. Most assumptions were obtained 
through personal communication with Yungelson. The binary 
fraction is set to 100%, while an IMF with exponent -2.5 is as- 
sumed. For the initial mass ratio distribution a broken power law 
is used, which is flat for close binaries. The distribution of initial 
orbital separations is taken to be logarithmically flat. A CE for- 
malism using the a parameter with a = 1 is used, however with 
an important change in the description: 

(M n + M zi ) (Mh - M lt ) (M u M 2 i M n M 2 i\ 



This equation yields eventual orbital periods which are typi- 
cally a factor of four larger than those obtained with Eq. (6). The 
eventually given DTD is for an instantaneous starburst, normal- 
ized to a formation of 4.7 M of stars per year. 



4.2. Han & Podsiadlowski 

lHan & Podsiadlowskil ([2004) only consider the SD channel for 
the formation of SNe la. Within this scenario, only progenitors 
which consist of a WD accreting from an unevolved or slightly 
evolved non-degenerate star are considered, i.e. the WD+MS 
channel in our terminology. They employ the Eggleton stellar 
evolution code for a Z = 0.02 population. A critical mass trans- 
fer rate toward the WD is given, above which no accretion takes 
place and mass transfer becomes non-conservative: 

[f) ^ = 5.3-10-^ ("wd- 0.4). (11) 

In this equation, X is the hydrogen mass fraction and masses 
and rates are expressed in M and M /yr respectively. Mass loss 
from the system is tak en to carry away the specific angular mo- 
mentum of the WD. A lMiller & Scalol(ll979l) IMF, constant ini- 
tial mass ratio distribution and logarithmically flat initial orbital 
separation distribution are adopted, while limits on mass, mass 
ratio and orbital period distributions are given. About 50% of 
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systems are binaries with an orbital period below 100 years. For 
CE evolution, the a-scenario is used with different values for a. 
A DTD is given for an instantaneous starburst of 10 11 M . No 
mass stripping is taken into account. 

4.3. Hachisu etal. 



In the same paper which describes the SD channel. [Hachisu et al.l 
(2008) use this scenario to carry out a population study them- 
selves. They only consider the SD scenario, through both for- 
mation channels, for a starburst with an exponent -2.5 IMF, flat 
mass ratio distribution and logarithmically flat orbital period dis- 
tribution. Unfortunately, only relative rates are given and the 
time binning is quite crude. As a result, only the morphologi- 
cal shape of this DTD can be compared to ours. 

Therefore, we also make a comparis on with the results of 
lHachisu et alJ dl999al) as presented by lHan & Pods iadlowski 
d2004l) . These do however not include the WD+RG channel 
(through which we do not obtain a significant number of SNe 
la anyway), nor the mass stripping effect for the WD+MS chan- 
nel. In this study the mass transfer rate toward the WD is taken 
to be 



dM 2 



At 



(12) 



In this equation, Tkh is the Kelvin-Helmholtz timescale, 
£rl the mass-radius exponent of the inner critical Roche lobe, 
and £ms that of the MS sta r. The resulting DTD is given by 
Han & Podsiadlowski (2004), using the same initial conditions 
as for their own model. 



4.4. Ruiteretal. 

The final study wit h which our results will be compared is that of 
iRuiter et all (120091) . This work investi gated both the SD a nd DP 
scenario using the evolution code of Belczvnski et al. (2008). 
This is done fo r a m etallicity of 0.02, a binary fraction of 50%, 
a lKroupa et al.l (1 19931) broken power law for the IMF with given 
mass range, a flat mass ratio distribution and an lAbtl (Il983l) ini- 
tial orbital separation distribution. For non-degenerate accretors, 
half of the mass that is transferred from the donor is accreted, 
while the other half is lost from the system with the specific an- 
gular momentum of the accretor. In the case of WD accretors, 
mass transfer is by default conservative, however the accretion 
rate is Eddington-limited. The ff-scenario for CE evolution is 
used, with discrete values for aA. The DTD is given both for an 
instantaneous burst of star formation and a constant star forma- 
tion rate during 10 Gyr. 

4.5. Discussion 

Adopting the elements of the respective codes described above, 
we carried out DTD calculations with the Brussels code. We find 
a good agreement with all of them, both in DTD shape and abso- 
lute number of SNe la. A comparison with the output obtained 
using our own assumptions (with (3 — 1 and the a-formalism) is 
shown in Figure 6 for the DD scenario and in Figure 7 for the 
SD scenario. There are a few differences, which are discussed 
below: 

• Our SD DTD shows a very high but short spike at very short 
de lay times, which in sh ape matches that of the WD+MS events 
by lHachisu et al. (2008). We do however not reproduce the sig- 
nificant number of events they find through the WD+RG chan- 



nel, since we do not have a non-negligible number of WD sys- 
tems traversing this progenitor region in the parameter space. 

• The aforementioned early spi ke is not present in the st udies 
of lHan & Podsiadlowskil d2004l) and lHachisu etal] dl999ah . As 
explained before, this is the result of the inclusion of mass strip- 
ping in our WD+MS channel, which was not done in those two 
studies. When we take c\ = 0, i.e. turn the effect off as well, 
the spike disappears and the early SD DTD is in good agreement 
with theirs. We do however retain a low background of SD SNe 
la beyond 1.4 Gyr which is not present in the comparison ones. 

• This low conti nuing SD rate is however confirmed by 
Ruiteretal. (2009), but these are the only events they find: the 
early maximum is not present is their s tudy. A possible explana- 
tion for this is that IRuiter et al.l ([2009 ) report that most of their 
SD events are created through the WD+RG channel (including 
subgiants) as opposed to the WD+MS channel, whereas in our 
case the opposite is true. It thus seems that the early maximum 
which is for some reason absent in their study corresponds to 
these events. 

• As far as the DD DTD is concerned, we f ind a very good 
agreem ent with both lYungelson & Livid d2000t) and IRuiter et al.l 
(2009). For the latter, this is however only true for fi = 0.5 
(which corresponds to their mass loss mechanism for RLOF to- 
ward a non-degenerate accretor) if we adopt a very small value 
for T] (< 0.25, simulating their assumption of specific accretor 
angular momentum loss). 

5. Further considerations and outlook 

In Sect. 3, a discrepancy in absolute number between theoreti- 
cal and observed SN la rates was discussed. A possible solution 
to this is stellar rotatio n: there seem to be indications (see e.g. 
lHabets & Zwaar] [l989) that it is common for components in bi- 
naries to rotate faster at birth than average single stars, which 
makes it possible that a significant fraction of primaries rotate 
faster than synchrono us. If so, their MS convec t ive core masses 
will b e larger (see e.g. Maeder & Mevnet 2000; Decressin et al.l 
2009). Hence, the same will be true for their remnant masses af- 
ter RLOF and the number of double WD binaries which in total 
mass exceed the Chandrasekhar mass will increase, leading to 
more SNe la through the DD scenario. We find that a convective 
core mass increase of only 10% is enough for the theoretical rate 
at 1 1 Gyr to match the observed one. The thus obtained DTD, for 
both the SD and DD scenario combined, is rendered in Figure 8. 
This figure also includes the v ery recent observational DTD as 
reported by Ma oz et al.l (l2010h . 

Also recently, Wang et al. (2010) proposed a SD mechanism 
in which instabilities in an accretion disk formed around the 
WD can lead to SNe la with very long delay times, of up to the 
Hubble time. However, the absolute number of events thus ob- 
tained at long delay times (beyond 3 Gyr) lies about an order of 
magnitude or more below that obtained through the DD scenario 
in the current paper, and hence even further below the observa- 
tions. As a result such mechanism, while beyond the scope of 
this work, does not alter our conclusions. At the other end of the 
DTD, the earliest SD events obtained in this paper are a result 
of the mass stripping effect and are therefore dependent on its 
efficiency. Another way to obtain such early SD even ts is a WD 
+ He star channel, as proposed bv iWang et al.l d20 09). However, 
even without including this channel, our computations are as- 
sured of very early events (at 50 Myr) through the DD scenario. 

The argument can also be extended to spiral galaxies, e.g. 
the Milky Way Galaxy. For the latter, it has been demon- 
strated in the past that the DD scenario can explain the Galactic 
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Fig. 6. DTD of a 10 6 M starburst obtained for the DD scenario with (3-1 and using the a-formalism for CE-evolutio n, compared 
with v arious other theoretical stud ies (see text). The solid l ine is the DTD from the present study, the dotted one is from lRuiter et al.l 
(2009) and the dashed one from Yungelson & Livid (120001) . 
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Fig. 7. DTD of a 10 6 M starburst obtained for the SD scenario with /3 — 1 and using the a-formalism for CE-evolutio n, compared 
with v arious other theoretica l studies (see text). Th e solid line is the DTD from the present study, the dotted one i s from lRuiter et al.l 
(2009), the dashed one from Hac hisu et alJ(ll999al) and the dashed dotted one from lHan &~P odsiadlowski (20041). 



birthrates well (see e.g. lHan et alj[l99l lHarJll998l) . The same 
has been done in previous stud ies with the Brussels code (see 
iDe Ponder & Vanbeverenir2004l) . which is an additional indica- 
tion that the difference between theoretical and observed rates 
in ellipticals may at least be partially caused by uncertainties on 
the observations in such galaxies. 



6. Conclusions 

We conclude that the single degenerate scenario alone cannot 
reproduce the observed distribution of delay times of type la su- 
pernovae. The double degenerate delay time distribution, possi- 
bly combined with the single degenerate one, does agree with 
the morphological shape of the observed distribution. However, 
assuming the a-formalism for common envelope evolution, this 
is only true if Roche lobe overflow is treated to occur quasi- 
conservatively. In this case, most type la supernovae are formed 



through such a canonical Roche lobe overflow followed by a 
common envelope phase, not through two successive common 
envelope phases. The absolute number of events critically de- 
pends on a whole range of parameters and descriptions, a fact 
which may be used to gain further knowledge about them once 
more detailed observations of this absolute number become 
available. At this time however, the theoretical models are un- 
able to reproduce the absolute number of observed type la su- 
pernovae, which is underestimated by a factor of at least three 
in elliptical galaxies with an age of 1 1 Gyr. This discrepancy is 
similar to those found in other population synthesis studies, but 
less pronounced in magnitude. It could be resolved by including 
the effects of rotation of binary components on their evolution. 
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